IEEE, 7° encuentro de Energía, Potencia, Instrumentación y Medidas. 16 y 17 de Octubre del 2008, Montevideo-URUGUAY. PAG 162 


Modelo Simulink de lámpara de vapor de sodio 
de alta presión (Na HID) para dimerizado en 
baja frecuencia 


José Inda, Gonzalo Casaravilla, Gabriel Barbat, Fernando Chiaramello 


Resumen-- Descripción y resultados de un modelo Simulink 
de lámparas de vapor de sodio, desarrollado para simular la 
corriente consumida por estas lámparas al ser alimentadas con 
una tensión sinusoidal de 50/60Hz con armónicos de baja y alta 
frecuencia, como por ejemplo, una onda recortada por un 
dimmer. Para determinar los parámetros del modelo, se mide 
simplemente tensión y corriente sobre la lámpara. 
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I. INTRODUCCIÓN 


La corriente consumida por una lámpara de vapor de sodio 
de alta presión tiene una forma de onda particular que 
necesita ser caracterizada si se quiere estudiar su consumo de 
armónicos. Este comportamiento no se puede simular con 
componentes lineales básicos. Si bien existen modelos para 
simular el consumo característico de corriente de una lámpara 
de vapor de sodio de alta presión (o más en general, cualquier 
lámpara de descarga), la mayor parte de los modelos 
desarrollados requieren de medidas difíciles de realizar (por 
ej. presión interna de la lámpara) para determinar los 
parámetros de simulación posteriores. También se 
encontraron modelos cuyos parámetros se obtienen mediante 
medidas más sencillas (midiendo tensiones y corrientes), pero 
se aplicaban al caso de fuentes sinusoidales [1] [4] o para 
fuentes de alta frecuencia (para simular un balasto 
electrónico, por ejemplo) [3], pero no se encontraron modelos 
que simulen el comportamiento de una de estas lámparas 
cuando son alimentadas por tensiones con componentes 
armónicos de relativamente baja frecuencia. Al controlar el 
flujo de una lámpara con un dimmer de tiristores, la tensión 
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de alimentación de la lámpara se vuelve una onda sinusoidal 
recortada, con componentes de Fourier de baja y alta 
frecuencia. Si se quiere estudiar el consumo de armónicos de 
corriente de una lámpara de descarga, los modelos recién 
presentados no sirven. En este trabajo se presenta un modelo 
Simulink de lámpara de vapor de sodio de alta presión que 
simula su consumo de corriente cuando es alimentada con 
una onda de tensión sinusoidal recortada, cuyos parámetros se 
obtienen midiendo tensión y corriente sobre la lámpara, y que 
se comporta bien a diferentes niveles de alimentación, con los 
mismos parámetros. 


II. BREVE DESCRIPCIÓN DE UNA LÁMPARA DE VAPOR DE SODIO 
DE ALTA PRESIÓN — COMPORTAMIENTO EN CONDI-CIONES 
NORMALES 


A efectos de limitar la corriente de consumo de una 
lámpara de descarga cuando se produce el arco, se debe poner 
en serie con la alimentación, un balasto (inductancia), como 
se ve en la fig. 1. 





Lámpara 


Fig. 1. Circuito de lámpara en funcionamiento con su balasto. 


Cuando la lámpara está en funcionamiento con una fuente 
sinusoidal, la tensión sobre la lámpara y la corriente que 
circula por ella toman la forma que se ve en la fig. 2. 

Las lámparas de descarga no tienen por qué ser 
alimentadas tensión sinusoidal. Los balastos 
electrónicos alimentan a estas lámparas con una onda 
cuadrada de alta frecuencia. 

En particular, se puede usar un controlador de flujo 
(dimmer) a tiristores, disminuyendo la tensión eficaz 
suministrada al recortar la forma de onda, y por lo tanto, 
disminuyendo el flujo lumínico producido por la lámpara, al 
mismo tiempo que la potencia consumida. 


con una 
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Fig. 2. Corriente y tensión registradas sobre una lámpara de vapor de sodio de 
alta presión en funcionamiento normal a 50Hz. 


Cuando la lámpara es alimentada con un dimmer como se ve 
en la fig. 3, las curvas de tensión y corriente resultan las que 
se ven en la fig 4. El objetivo del modelo presentado aquí es 
que tenga un comportamiento como el descrito. 





Fig. 3. Lámpara de vapor de sodio de alta presión alimentada con un dimmer. 
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Fig. 4. Tensión y correiente registradas sobre la lámpara alimentada con un 
dimmer 


HI. MODELO 


Una versión previa al modelo propuesto en este trabajo es 
el mostrado la fig. 5. 

A la izquierda se encuentra la parte del circuito eléctrico y 
consta de un amperímetro y una fuente controlable de tensión. 
La parte de señal del circuito, que se encuentra a la derecha, 
toma la medida del amperímetro y calcula la tensión que 
impondrá a la fuente controlable. 


Lr] 


Amperímetro 






Fuente controlable 
de tensión 





Conn2 
Fig. 5. Circuito de modelo previo 


En detalle, la parte de señal calcula una resistencia en 
función de la corriente 


v 
R(¿)=min Rea 0) 


donde Rsar y Vc son parámetros del modelo de la lámpara. 
La señal R(i) pasa por un filtro de transferencia 


1 
H(s)+—=— (2) 


(a. s +1)? 


para obtener así la “resistencia instantánea” de la lámpara. 

Esta señal es multiplicada por la corriente para obtener la 
señal de tensión que se ingresa a la fuente controlable de 
tensión. 

Este modelo tiene 3 parámetros 
comportamiento: Vc, Rsat Y a. 

Rsa y a son fijos con lámpara mientras que Vc depende de 
la potencia media consumida por la lámpara. A tales efectos 
se agregó al modelo un bloque que ajusta dicha dependencia. 
El modelo final propuesto se muestra en la fig. 6. 

El bloque que se agregó calcula 


que afectan su 


P=v+i (3) 


De esta señal se integra el último ciclo (los ultimos 
20ms/16.7ms) para hallar la potencia media consumida P,, y 
se calcula Vc a partir de P,, como 


=A-P +B 
Ve ia: 4) 


donde A y B son parámetros del modelo de la lámpara. 
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Fig. 6. Circuito del modelo propuesto. 


IV. JUSTIFICACIÓN DEL MODELO 


En primer lugar, este modelo calcula v en función de i (y 
no a la inversa) porque en serie está el balasto, que impone la 
corriente sobre la lámpara más que su tensión. Esto es lo 
usual en otros modelos vistos de lámpara. 

En segundo lugar, el paso intermedio de calcular R(i) y 
pasarlo por una transferencia H(s) de 2 polos fue lo que mejor 
funcionó de todo lo que se probó para reproducir el sobretiro 
de tensión característico asociado con el transitorio 
establecimiento del arco eléctrico. De todos modos, algunas 
características favorecieron la elección de esta alternativa. 

Si se observan las curvas de tensión y corriente sobre las 
lámparas, se ve que ambas se anulan al mismo tiempo, lo que 
representa una carga completamente resistiva (en contraste 
con cargas  inductivas O  capacitivas) aunque no 
necesariamente sea una resistencia constante en el tiempo. 

Además, se estudió la resistencia (variable) de la lámpara 
en función del tiempo, mostrando que es razonable elegir 
calcular la resistencia instantánea como paso intermedio, y 
como resultado de una transferencia que imponga algún 
retardo. 

Otras alternativas fueron hacer la resistencia en función de 
la potencia instantánea R(P) y pasarlo por una transferencia, 
similar a como estaba planteado en [3], o calcular la tensión 
directamente a partir de la corriente v(i) y pasarlo por una 
transferencia. También se probaron otras transferencias y se 
encontró que la transferencia de la ecuación (2) dio los 
mejores resultados. 

Para determinar la función R(i) se tomó en cuenta que la 
forma de onda de la tensión es similar a una onda cuadrada, 
manteniéndose casi constante durante algo menos de medio 
período, en los momentos en que la corriente toma valores 
absolutos más grandes. Por lo tanto, se puede decir que para 
corrientes relativamente grandes se cumple 
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donde Vc es la tensión sobre la lámpara cuando ésta se 
mantiene constante. 

Sin embargo, cuando la corriente es pequeña, no se cumple 
(5), haciéndose R más pequeña. 

De esta manera se determinó la ecuación (1). El parámetro 
Vc corresponde a tensión mencionada en (5), y Ra, es un 
valor tope de resistencia en principio arbitrario. 

R(i) se pasa entonces por la transferencia que la suaviza y 
la retrasa produciendo al final la subida que tiene la forma de 
onda de la tensión cada vez que cambia de signo. 

Al comparar el modelo con las diferentes medidas se 
dedujo que los parámetros Ry, y a solo dependen de la 
lámpara y que Vc también depende de las condiciones de 
alimentación. 

Al comparar las formas de onda de tensión sobre la 
lámpara en condiciones en que la potencia media consumida 
Pm fuera la misma, se observó que Vc (la parte plana de las 
mismas) se mantiene similar, y si Pm cambia, Vc cambia 
sensiblemente. 

Finalmente, al observar la dependencia de Vç con Pm se 
encontró que se trata de una dependencia lineal del tipo de la 
ecuación (4). 

Comparando las medidas con los resultados de la 
simulación se ajustaron todos los parámetros y los resultados 
fueron satisfactorios. 


V. RESULTADOS 


Se hicieron un conjunto de medidas aparte de las que se 
usaron para determinar el modelo y en todos los casos los 
resultados fueron satisfactorios. En particular se agregó un 
capacitor de compensación de reactiva entre la fuente (que 
tiene cierta impedancia de cortocircuito) y el dimmer como es 
mostrado en la fig. 7, y se midieron las formas de onda de 
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tensión y corriente de línea y sobre la lámpara. 










Dimmer Balasto 


Lámpara 


Fig. 7. Circuito con capacitor de compensación de reactiva. 


El comportamiento del modelo es mejor en corriente que 
en tensión, que por otro lado, es donde usualmente se necesita 
que sea mejor. En la fig. 8 se observa que prácticamente no se 
perciben diferencias entre la forma de onda de corriente 
medida y obtenida en la simulación, cuando la fuente es 
sinusoidal. 

En las fig. 10 y 11 se observan las formas de corriente y su 
contenido armónico, cuando la alimentación es una onda 
recortada. 
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Fig. 8. Corriente sobre la lámpara, medida y simulada con el modelo, para el 
caso de fuente sinusoidal. 
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Fig. 9. Tensión sobre la lámpara, medida y simulada con el modelo, para el caso 
de fuente sinusoidal. 
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Fig. 10. Corriente sobre la lámpara, medida y simulada con el modelo, para el 
caso de fuente sinusoidal recortada. 
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Fig. 11. Armónicos de corriente sobre la lámpara, para el caso de fuente 
recortada 


Cuando se agrega el capacitor de compensación de 
reactiva, se observa una gran distorsión armónica en la red, 
debido a que se produce resonancia entre dicho capacitor, y la 
inductancia de cortocircuito de la red, lo cual es reportado 
(sin dimmer) en [2]. En el anexo se hacen las cuentas que 
muestran la dependencia de la frecuencia de resonancia en 
función de la relación entre la potencia aparente nominal del 
transformador de distribución de la red, y la potencia aparente 
de la instalación de lámparas, obteniéndose para este caso una 
frecuencia de resonancia de 255Hz, muy cercana el armónico 
5. Al simular este efecto con el modelo de lámpara presentado 
aquí, se obtienen también buenos resultados como se ve en las 
fig. 12, 13 y 14. 
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Fig. 12. Corriente sobre la red, medida y simulada con el modelo, para el caso 
con capacitor de compensación de reactiva. 
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Fig. 13. Tensión de alimentación del conjunto dimmer-balasto-lámpara, medida 
y simulada con el modelo, para el caso con capacitor de compensación de 
reactiva. 
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Fig. 14. Armónicos de corriente sobre la red, medida y simulada con el modelo, 
para el caso con compensación de reactiva. 




















VI. CONCLUSIÓN 


Se determinó un modelo Simulink de lámpara de vapor de 
sodio de alta presión, con comportamiento satisfactorio 
cuando la tensión de alimentación es no puramente sinusoidal 
y tiene armónicos de baja frecuencia. 

Si bien existen varios modelos de este tipo de lámparas 
para simulación, la mayoría o bien no cubren el caso 


particular de alimentación recién mencionado, o bien 
requieren medidas complicadas. 

El modelo aquí presentado cubre satisfactoriamente el caso 
recién mencionado, además de tener lo posibilidad de hallar 
sus parámetros midiendo solamente tensiones y corrientes 


sobre la lámpara. 


VII. APÉNDICE 


En una instalación de lámparas de descarga con 
capacitores de compensación de reactiva se produce un 
fenómeno de resonancia que amplifica respuestas no lineales 
de las lámparas. La frecuencia de resonancia depende de los 
capacitores de compensación de reactiva y la inductancia de 
cortocircuito del transformador de distribución aguas arriba, 
que dependen a su vez, de la potencia de la instalación de las 
lámparas y la potencia nominal del transformador. En este 
apéndice se encuentra una expresión para dicha dependencia 
en una instalación típica, y se muestra que la frecuencia de 
resonancia es menor (y por lo tanto peor) cuando la potencia 
nominal del transformador es cercana a la potencia aparente 
de la instalación de lámparas, llegando a ser la frecuencia de 
resonancia, comparable a la frecuencia fundamental de la red. 

La frecuencia de resonancia está dada por la ecuación 


A 1 1 (6) 
r 2.1 LC 


donde L es la impedancia de cortocircuito del transformador 
de distribución, y C es la capacidad de fase de la instalación. 

Considerando la componente fundamental de la carga, 
dada una determinada instalación equilibrada de lámparas, se 
conoce su potencia de consumo P, y su potencia aparente S (y 
por lo tanto su factor de potencia cos). La potencia P 
depende únicamente de las lámparas, mientras que S también 
depende de los balastos y los capacitores de compensación de 
reactiva. 

En general, dada una instalación con solamente las 
lámparas y los balastos, se conoce, al igual que su potencia P, 
su potencia aparente Sz (y por lo tanto su factor de potencia 
cos QL), y se eligen los capacitores para obtener una S deseada 
en el punto de facturación. 

La capacidad de fase de la instalación resulta determinada 
por: 


C= (7) 
gü” 


donde O, es el consumo de potencia reactiva del conjunto 
balasto lámpara sin compensar, Q es el consumo de potencia 
reactiva deseada luego de compensar, w es la frecuencia 
angular de la red y U es la tensión nominal de la red (entre 
fases). 

Operando se obtiene 
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o=15? - P? (8) 
0, =5;-P? (9) 
P = S -coso (10) 


De esta forma se obtiene C en función de S: 


S 
C- Sopr 1- cos? p, Sd (11) 


oU? cos py 


Por otro lado, se puede expresar la inductancia de 
cortocircuito del transformador de distribución en función de 
su potencia nominal Sz: 











2 
U 
CL X cC (12) 
o: Sy 
donde Xcc es la impedancia de cortocircuito en p.u. 
Finalmente, sustituyendo en la eq. (6) se obtiene 
FOIS 
E as 
fo CC <52 fi- eos? m ia 
cos py L 


que muestra la relación entre la frecuencia de resonancia de 
la instalación, la potencia aparente de la misma, y la potencia 
nominal del transformador de distribución aguas arriba. 

En el ensayo de laboratorio, montado con luminarias 
comunes, el factor de potencia sin compensar resultó 





cosp, =0,45 (14) 
Si el factor de potencia compensado deseado es 
cos p=0,92 (15) 
sustituyendo en (13) se obtiene 
i =0,84 - ‘r . = (16) 
0 CC 


Evaluando este resultado, suponiendo Xcc en valores 
normales de 5% y 10% se observan las curvas de la fig. 15. 
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Fig. 15. Frecuencia de resonancia de una instalación de lámparas en función de 
la potencia aparente consumida y la potencia del transformador de distribución, 
para dos valores de impedancia de cortocircuito del mismo. 


En el ensayo del laboratorio, los valores obtenidos fueron: 


=0,4 
cos o L 0,45 (17) 
cos p=0,85 (18) 
S = 2049 - VA (19) 
Sy =2500 - VA (20) 
X cc =5% Q1) 
-E=5,1 (22) 

fo 
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